Modulation structure of a frustrated spin ladder 
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We study a two- leg spin- 1/2 ladder with isotropic exchanges and biquadratic interactions in the 
basic plaquettes. It is shown that for the extremely frustrated case, the system exhibits a self- 
organized phase separation. In some parameter regions, the singlet rungs form a Wigner-like lattice 
in the triplet-rung host. There are three types of elementary excitations in this modulation phase, 
i.e., the spinous in a triplet domain, the broken singlet rungs and the deformation of the Wigner-like 
lattice. The flux phase induced by an external magnetic field in the rung-dimerized phase is also 
discussed. 
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After the discovery of the ladder superconductors, the current interest in the coupled spin chains has been greatly 
renewed.lj It is well established that the regular spin-1/2 ladders with even number of legs have a gaped spin-liquid 
ground state, while the odd-legged ladders have a gapless spin-liquid ground state. On the other hand, the generalized 
spin ladders including other couplings beyond the nearest neighbor exchanges, which can interpolate among a variety 
of systems, have been demonstrated exhibiting remarkably rich behaviorE~l3 The hamiltonian of the generalized spin 
ladder reads: 

1 ^ 1 ^ 

1 ^ 1 ^ 

1 ^ 1 ^ 

where aj {tj) are Pauli matrices on site j of the upper (lower) leg, Jq {a — 1, 2, 3) are the exchange coupling constants 
and Ua are the biquadratic coupling constants. The hamiltonian (1) contains all possible S't/(2)-invariant interactions 
in a basic pilaquette formed by two nearest rungs. We note the biquadratic interactions may be, effectively mediated 
by phononsEI and their importance for some properties of CuO^ plaquette has been pointed out.l3 In addition, current 
experiments have revealed that such multi-spin interactions are realized in the two-dimensional (2D) solid ^He)^ 2D 
Wigner solid of electrons formed in a Si inversion layerO, and the hcc solid ^He^ 

The model (Ijuhas been studied bv^aiany authors in some parameter regions. The possible phases including the 
dimerized phase J3 the Haldane phase aB and the rung-dimerized phase have been reported. In addition, the model 
(1) has a variety of integrable caseal3'EJ and some gapless phases have been found. All these results indicate that the 
generalized spin ladder has a very rich ground state phase diagram. In this letter, we study the extremely frustrated, 
i.e., Ji = J2, Ui = U2 case. We show that there are three possible phases in this case, i.e., a rung-dimerized (RD) 
phase, a triplet-rung (TR) phase interpolating between the gapless spin liquid and the Haldane spin liquidEJ or the 
VBS statejlj and a new gapful phase with modulation structure. The latter phase is a mixed state (MS) consisting of 
both singlet rungs and triplet rungs and its structure is very similar to a Wigner lattice (crystallization of the singlet 
rungs) or the flux phase of the type II superconductors (with the singlet rungs as the effective fluxes). 

In the extremely frustrated case (Ji — J2, Ui — U2), the hamiltonian (1) can be rewritten as 
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N N 

+C/^(2 - S]){2 - 5|+i) -JY.{2- Sf) + C, (2) 

where Sj = {dj+T.j)/2 is the total spin of the j-th rung, J 2Ji + C/i, (3 = -2C/i/(2Ji + C/i), U = U3, J = J3 + U3-2U1 
and C = J3/2 + J73/4 — 5C/1/2 is an irrelevant constant. For convenience, we put J = 1 and omit the irrelevant constant 



1 



C in the following text. The total spin Sj takes two possible values 0, 1, corresponding to the singlet rungs and triplet 
rungs, respectively. Based on Eq.(2), the phase diagram of the system can be conjectured roughly. There are three 
possible types of ground state configurations. For a large enough J, the singlet rungs are more stable than the triplet 
rungs. The ground state is a simple product of TV singlet rungs, i.e., a rung-dimerized state. For a small J and large 
U , the triplet rungs are dominant over the singlet rungs and the ground state is exactly the same of a generalized 
spin-1 chain. In some intermediate J, U regions, the singlet rungs and the triplet rungs may coexist in the ground 
state. Since there is no genuine interactions between the singlet rungs and the triplet rungs as we can read off from 
Eq.(2), a self-organized phase separation occurs in this case. The singlet rungs simply cut the system into disconnected 
domains of triplet rungs and a positive U will stablize a Wigner-like lattice of the singlet rungs in the triplet liquid. 
Each triplet domain in the mixed state behaves as an open spin-1 chain. Such a mixed state can only be accompanied 
by a negative boundary energy, as the flux phase in a type II superconductor. We note that [3r-^ ±1 represent two 
integrable points of the model (2) and when /3 = — 1/3, a VBS ground state can be constructed.t2l Very interestingly, 
the translational invariance is broken in the MS phase and the ground state is no longer a spin liquid but a periodic 
array of spin-liquid domains. 

To give a clear picture, we study one of the integrable cases (3 — —1. In this case, each triplet cL©iM|in behaves as an 
S'J7(3)-invariant spin-1 chain with open boundaries and can be solved exactly via Bethe ansatzJlil't^ Let us consider 
a triplet domain with M rungs. The effective hamiltonian of this domain reads 
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with Sj — 1. The solution of Eq.(3) is exactly the same of the S'[/(3)-invariant t — J modelEl with pure open 
boundaries. Its spectrum is determined by the following Bethe ansatz equations (BAE's) 
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where Xj and are the rapidities of the spinous and M2 < Mi < M. The eigen energy of the triplet domain reads: 
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Notice that the second term in Eq.(5) is also relevant since the number of triplet rungs is not fixed in the whole 
system. It contributes an amount of -2 to the boundary energy, which is crucial to stablize an MS. For M >> 1, the 
ground state energy of the triplet domain can be expressed as Em = Mcq + eb{M), where Cq is the energy density 
of an infinite SU{3) spin-1 chain, eb(M) is the boundary energy including the 0{M~^) finite size correction. Denote 
the ground state distributions of A and /i for Af — > 00 as pi{X) and P2(m) respectively. From the BAE's (4) we get 
the following integral equations 



Pi(A) + [2]pi(A)-[l]p2(A) 
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a2(A) -I- ai(A), 



P2(m) + [2]P2(m) - [1]Pi(m) = ^"2(m), 
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where a„(A) = n/27r[A^ -I- (n/2)^] and «J is, the integrable operator with kernel a„(A). By solving the above integral 
equations with Fourier transformation|lJil2l we obtain 
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where p°(A) is the density of A of the ground state for M ^ oo with periodic boundary conditions. The 5IX\ term in 
the third equation of Eq.(7) comes from the A = mode, which is forbidden in an open boundary system.tSEJ cq and 
eb(M) can be derived as 

eo = _ / /i(^dA + 2 = 2-ln3 " 



A2 + 1/4 3V3' 
^b{M) = - / 2 + O(M-i) = ^^/37^ - 4 + 0{M-^). (8) 

Based on Eq.(8), the phase boundaries can be determined exactly. As we discussed above, there are three possible 
phases, i.e., the RD phase, the TR liquid and the MS phase containing both singlet rungs and triplet rungs. We note 
that the density of the ground state energy of the rung-dimerized state reads = 4[/ — 2 J, as we can easily derive 
from Eq.(2). Therefore, the phase boundary between the RD phase and the TR liquid phase is given by 4U — 2J = to- 
On the other hand, generating a singlet rung in the TR liquid implys a broken triplet rung and an open boundary. 
The excitation energy of this process is eb(oo) — eo — 2 J. Therefore, the phase boundary between the TR phase and 
the mixed state is given by 2J = eb(oo) — eg. We note the singlet rungs are very similar to the fluxes in a type II 
superconductor. Here eb(oo) < corresponds to the surface energy, while — 2J — eo corresponds to the self energy of 
the fluxes. When 2 J > eb(oo) — eg, some singlet rungs appear in the triplet liquid. Here Jd = (eh(cx)) — eo)/2 serves 
as the lower critical field. Due to the finitepSiza. correction of eb{M), the singlet rungs are unfavorable to close each 
other (note ef,(M) is an increasing functior£3il3 of 1/M). The lengths of the triplet domains are mainly controlled 
by the finite size correction. Suppose we have Nm — N/{M + 1) triplet domains with lengths AI + 5m respectively. 
The small shifts 5m satisfy the condition 5m = 0. The total correction energy reads 

Nm Nm , 

E, = Y, + 5 m) - e6(oo)] ^ j2 jTTs::,- 

By minimizing Eq.(9), we get 5m = 0, which indicates a periodic array of the singlet rungs in the whole system. 
This novel modulation structure is very similar to a Wigner lattice or an Abrikosov lattice but with a very different 
physical interpretation. As in the type II superconductors, there is also an upper bound of J corresponding to the 
upper critical field, which determines the phase boundary between the RD phase and the MS phase. Suppose an MS 
is stable at M = M{J). The phase boundary is given by 

4(7 - 2J - eM(j) + Jl^Jf^l^MiJ) - 2 J]. (10) 

When M{J) ^ 1, J^2 ~^ 41/. 

Now we turn to the elementary excitations in the MS phase. There are three types of elementary excitations in 
this case, i.e., the spinous in the triplet domains, broken singlet rungs and the deformation of the Wigner- like lattice. 
Suppose the lattice is stablized with a period M + 1 {M triplet rungs and one singlet rung). The spinous in the present 
case have a finite energy gap in the order of due to the finite length. oLthe domains. For M » 1, the energy 
gap can be derived with the well known finite-size-correction techniques.EZtiliAs in the usual integrable models, the 
excitation energies of the spinons are given by the so-called dressed energiest2l which in our case read (for M — > oo) 

ei(A) = -2™i(A) + [l]e2(A) - [2]ei(A), 

e2(Ai) = Wl{^i) - [2]e2(M)- (11) 

The velocities of the two branches of low-energy spinons are equivalent, which read 

y = lim ^ lini (12) 

The energy gap associated with the spinon excitations is thus A w t:v/M. There are three different pictures of the 
spinon excitations: (i)For M = 3n, there is no holes in the A- and /i-seas in the ground state. The spinons are 
generated by spin flips. The simplest excitation is a one A-hole and two /Lt-hole state. (ii)For M ^ 3ri, there are some 
A- and ^- holes in the ground state. The mobility of the holes gives the simplest excitation. (iii)For arbitrary M, 
there are string excitations. We note all these excitations have the same energy gap. 

A broken singlet rung in the modulation phase indicate that two neighboring triplet domains combined to a 
single domain with length 2M + 1. In a stable MS phase, M is determined by minimizing the per site energy 
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CM = [eb(M) - eo - 2J]/(M + 1) + eo- For large M, from dcM/dM = we readily obtain J - [e6(oo) - eo]/2 - Af-^. 
Therefore, the excitation energy of a broken singlet rung is eh(2M + 1) — eh(M) + eo + 2J, which is still in the order 
of M-\ 

The deformations of the Wigner-like lattice represent another type of excitations in the modulation phase. This 
type of excitations are static rather than dynamic due to the complete phase separation. Suppose one singlet rung 
is moved from its equilibrium position by SM . The two neighboring triplet domains connected by this singlet rung 
are thus enlarged and compressed by 6M respectively. The energy of this process comes mainly from the finite size 
corrections eb{M + SM) + eb{M - SM) - 2et{M) and is in the order of M'^ (notice that eb{M) - eh(oo) ~ M"!), 
implying the larger the M, the more unstable the lattice. 

A ladder system may exhibit interesting behavior in an external magnetic field. In the TR phase, the system 
behaves as a two-component Luttinger liquid. The response of the system to the magnetic field is rather usual, i.e., 
the zero field susceptibility shows a simple Pauli law. In the MS phase, a magnetic field will enlarge the period of the 
Wigner-like lattice and at a critical field He, the singlet rungs are no longer stable even in the ground state, implying 
a phase transition between TR phase and MS phase. Two different situations may appear when a magnetic field is 
applied on the RD phase. Roughly speaking, the magnetic field depresses the effect of J. Therefore, with the increase 
of the field, the system flows either toward the TR phase or toward the MS phase. If the system flows to the MS 
under a magnetic field, some flux phase may appear. Here the triplet rungs with Sz — ^ serve as the fluxes in the 
singlet- rung host. Since the "repulsive interaction" U occurs only between the nearest neighbor singlet rungs, the only 
possible flux phase has the structure |0 >i >2 (^\0 >3 (8) • • • (E" |0 >n-i (^\^ >n, where |0 >i and |1 >j indicate the 
singlet rungs and the triplet rungs, respectively and N (even) represents the length of the ladder. With the increase 
of the external field, the fluxes will clusterize and finally form a triplet-rung liquid or a completely polarized state. 

In conclusion, a generalized spin ladder is studied. It is found that in the extremely frustrated case, a modulation 
structure which represents a Wigner-like lattice of the singlet rungs in the triplet-rung host can exist in some parameter 
regions. Though only the /3 = — 1 case is studied in detail, similar phenomena may exist for arbitrary (3. The only 
difference is that the TR domains in /3 = — 1 case may be replaced by Haldane domains or VBS domains. 
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